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Abstract 

We consider pair production of new down- type isosinglet quarks originating from Ea, which is the favorite gauge 
symmetry group in superstring inspired GUT models. The study concentrates on the possibility of observing the pair 
production of the lightest of the new quarks, D, in the ATLAS detector at the forthcoming LHC accelerator, in the 
channel DD — > ZjZj . Both signal and background events are studied using tree level event generators based on Monte 
Carlo techniques. The detector effects are taken into account using the ATLAS fast simulation tool, ATLFAST. It is 
shown that ATLAS can observe the D quark within the first year of low luminosity LHC operation if its mass is less than 
650 GeV. For the case of two years of full luminosity running, 1 TeV can be reached with about three sigma significance. 



1 Introduction 

If observed, the long awaited discovery of the Higgs particle at the LHC experiments dEI will complete the validation of 
the basic principles of the Standard Model (SM). However, the well known deficiencies of the SM, such as the arbitrariness 
of the fermion mass spectrum and mixings, the number of families, the real unification of the fundamental interactions 
and the origin of baryon asymmetry of the universe require extensions of the SM to achieve more complete theories. In 
general, these extensions predict the existence of new fundamental particles and interactions. The forthcoming LHC will 
give the opportunity to explore new colored particles and their interactions to test these predictions. Three types of new 
quarks (see |;3| for a general classification) are of special interest: the fourth SM family quarks, up type and down type 
weak isosinglet quarks. The existence of the fourth SM family is favored by flavor democracy (see |4l and references 
therein for details), Q=2/3 quarks are predicted by the little Higgs model 1 5 1 and Q = —1/3 isosinglet quarks are predicted 
by Grand Unification Theories (GUTs), with Eq as the unification group |6 |. The GUT models permit solving at least two 
of the above mentioned problems, namely, the complete unification of the fundamental interactions (except gravity) and 
the baryon asymmetry of the observed universe by merging strong and electroweak interactions in a single gauge group. 
Theories adding gravity to the unification of fundamental forces, superstring and supergravity theories |7 | also favor Eq 
as a gauge symmetry group when compactified from 10 (or 11) dimensions down to the 3+1 that we observe (see O and 
references therein for a review of Ee GUTs). 

For LHC, the production and observation of the first and second type of new quarks have been investigated in |Ij i9j and 
ISl llOl respectively. In this work, we study the possibility to observe the third, down type isosinglet quarks predicted by 
the Eq-GVT model, at LHC in general, and, specifically in the ATLAS experiment] 1 1 using the 4 lepton and 2 jet channel. 
The current experimental limit on the mass of an isosinglet quark is m > 199 GeV il U . The detailed study for down type 
isosinglet quark signatures at the Tevatron has been recently performed in 1121 and it has been shown that the upgraded 
Tevatron would allow a mass reach up to about 300 GeV. 



2 The Model 



If the group structure of the SM, SUc'{3) x SUw{2) x Uy{1), originates from the breaking of the Eq GUT scale down 
to the electroweak scale, then the quark sector of the SM is extended in the following manner; 
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As shown, each SM family is extended by the addition of an isosinglet quark. The new quarks are denoted by letters D, 
S, and B. The mixings between these and SM down type quarks is responsible for the decays of the new quarks. 

These mixings increase the number of angles (A^e) and phases (A^$) with respect to the SM CKM matrix. For a general 
multi-quark model, one has 1131 : 

Ne - Nx{l + m-^^^^) (2) 
= {N-l)x{l + m- ^ ) 

where I and m are the numbers of the up-type and down-type quarks respectively; and N is the number of SU{2)w 
doublets formed by left handed quarks. In the case of the Eq model, we have m = 21 = 2N = 6 and Eq. (|2ji yields 
A^e = 12 and iV$ = 7 . The special case m = Z + l = iV+l = 4, considered in 1121 . yields A^e = 6 and iV$ = 3 which 
coincides with the number of parameters in the Little Higgs models with one additional isosinglet up-type quark |j5j| . 

In this study, the intrafamily mixings of the new quarks are assumed to be dominant with respect to their inter-family 
mixings. In addition, as for the SM hierarchy, the D quark is taken to be the lightest one. The usual CKM mixings, 
represented by superscript 9, are taken to be in the up sector for simplicity of calculation (which does not affect the 
results). Therefore, the Lagrangian relevant for the decay of the D quark becomes L14J ; 



Cj, = ^^J^^[MS„(l-75)dcos0 + Aa(l-75)i?sin0]M^« (3) 
2\/2smOw 



4sin( 



'w 



sm0cos0- 

7 d-ja (1 - 75) D 

cos Ow 



+ h.c. 



The measured values of Vud,Vus,Vub constrain the d and D mixing angle to | siii0| < 0.07 assuming the squared 
sum of row elements of the new 3x4 CKM matrix give unity (see tlTl and references therein for CKM matrix related 
measurements). The total decay width and the contribution by neutral and charged currents were already estimated in 1 14J . 
As reported in this work, the D quark decays through a W boson with a branching ratio of 67% and through a Z boson 
with a branching ratio of 33%. The total width of the D quark as a function of its mass is shown in Fig.^for the illustrative 
value of sin = 0.05. It is seen that the D quark has a rather narrow width and becomes even narrower with decreasing 
value of (j) since it scales through a sin^ (j) dependence. If the Higgs boson exists, in addition to these two modes, D quark 
might also decay via the D —* H d channel which is available due to Z? — d mixing. The branching ratio of this channel 
for the case of niH = 120 GeV and sin0 = 0.05 is calculated to be about 25%, reducing the branching ratios of the 
previously discussed neutral and charged channels to 50% and 25%, respectively I12II15I . However, this study will not 
take into account the possible existence of the Higgs boson and will concentrate on the pair production of the D quarks 
which is approximately independent of the value of sin (j). One should note that the additional consideration of single 
production of D quark increases the overall D quark production rate; for example by about 40% for rriD = 800 GeV if 
sin <j> = 0.05. However, final state particles and SM backgrounds are different from the pair production case, making it a 
rather different process to study. 



3 Pair Production at LHC - signal at generator level 

The main tree level Feynman diagrams for the pair production of D quarks at LHC are presented in Fig. |2l The gDD 
and jDD vertices are the same as their SM down quark counterparts. The modification to the ZDD vertex due Xo d — D 
mixing can be neglected due to the small value of sin (j). The Lagrangian in Eq. (|3} was implemented into tree level 
event generators, Comphep 1161 version 4.3 and Madgraph 1171 version 2.3 . The total pair production cross sections 
from these two Monte Carlo generators are shown in Fig. |3] The difference in cross section calculated by these two 
generators, the first one, based on full matrix element calculation and the other on the numerical methods is less than 5% 
for the range of D quark mass from 400 to 1400 GeV. The impact of uncertainties in parton distribution functions (PDFs) 
1181 . is calculated by using different PDF sets, to be less than 10% for the same range. For example at 771 = 800 GeV 
and ~ m\, the cross section values are 450 (CompHep, CTEQ6L1) and 468 (CompHep, CTEQ5L) versus 449 
(MadGraph, CTEQ6L1) and 459 (MadGraph, CTEQ5L) fb with an error of about one percent. For the same PDF set. 
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Figure 2: The tree level Feynman diagrams for the pair production of isosinglet quarks 
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Figure 3: The DD pair total production cross section (solid line) as a function of D quark mass is shown as calculated 
with CompHep (triangles) and with MadGraph (circles). The dashed line is for the gluon contribution and the dotted line 
is for the quark contribution. The two generators agree within 2% and the effect of PDF uncertainties is computed to be 
less than 10% through 5 orders of magnitude in the cross section. 
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the two programs give the same answer validating both of the MC generators against each other and the implementation 
of the model. The same figure also contains the partial {gg and qq) contributions showing that the largest contribution 
to the total cross section comes from the first three diagrams for D quark masses <1 100 GeV, while for higher D quark 
masses, contributions from s-channel qq subprocesses becomes dominant. For these computations, qq are assumed to be 
only from the first quark family since, the contribution to the total cross section from ss is about 10 times smaller and the 
contribution from cc and bb are about 100 times smaller. The <-channel diagrams mediated by Z and W bosons, shown 
on the bottom row of Fig.|2j which are suppressed by the small value of sin (f> (for example 0.4 fb at mj^ = 800 GeV) are 
also included in the signal generation. 

The isosinglet quarks being too heavy are expected to immediately decay into SM particles. The decay channels leading to 
possible discovery for the D quark pairs are summarized in Tabled We have initially focused on the 4 lepton final states 
of the neutral channel only: although it has the smallest branching ratio, the possibility of reconstructing the invariant 
mass for Z bosons and thus for both D quarks makes it favorable for a feasibility study. Therefore the final state we 
would be looking for is composed of two high transverse momenta jets and two Z bosons, all coming from the decay of 
the D quarks. The high transverse momentum of the jets coming from the D quark decays can be used to distinguish the 
signal events from the background. 

All the SM processes allowed in proton collisions yielding two Z bosons and two jets (originating from any parton except 
b and t quarks) were considered as background. These can be classified in three categories via the initial state partons: 
qq, qg and gg where g stands for gluon and q can be any quark or anti-quark from the first two families. The contribution 
from third quark family is assumed to be negligibly small due to mass and PDF arguments. Although a complete list 
can be seen in l!20';1, one should note that the 78% of the background cross section originates from the processes in first 
two categories : qg qgZZ (where q ~ u, d, d, c, s, s) contributes 58%, and qq ggZZ (where q = d, u, s, c) 
contributes 20%. The simple requirements imposed at the generator level are: 



\Vp\ < 2.5, 

Pt,p > 100 GeV 

Rpp > 0.4 , 

\vz\ < 5.0 



(4) 



where R is the cone separation angle between two partons [p =d,d), ijp and r/z are pseudorapidities of a parton and Z 
boson respectively, and Pr^p is the parton transverse momentum. The selection of the jy region is driven by partonic 
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Table 1: The promising signal channels. The fourth column contains the branching ratios of the SM particles, whereas 
the last column has the total branching ratios. 



DD 


Final State 


Expected Signal 


Decay B.R. 


Total B.R. 


Z Zdd 
0.33 X 0.33 


Z 11 z a 


41 +2 jet 


0.07 X 0.07 


0.0005 


Z ^IIZ ^ vv 


21 +2jet + 


2 X 0.07 X 0.2 


0.0028 


Z 11 Z qq 


21 +4 jet 


2 X 0.07 X 0.7 


0.0107 


ZWdu 
2 X 0.33 X 0.67 


Z IIW Iv 


31 +2jet + ^T 


0.07 X 0.21 


0.0065 


Z 11 W qq 


21 +4 jet 


0.07 X 0.68 


0.0211 



spectra pseudorapitidy distributions, which are peaked in the barrel detector area. The signal cross section was calculated 
with both generators as a function of the D quark mass, but the background only with Madgraph as it is faster in numerical 
evaluation. For mo ~ 800 GeV, using the generator level cuts listed in Eq. (|4j, the cross section of the ZZ 2 jet channel 
is found to be (Tsignai = 45.4 ± 2 fb (CompHep) whereas the SM background for the same final state particles is 
(Jbg = 345 ± 17 fb (MadGraph). Already at this level, a significance of S/^/B w 2.4 can be obtained with an integrated 
luminosity of 1 fb~^. 



4 Observation in ATLAS detector using 4/ 2 jet channel 

Using the preselection cuts listed in the previous section, 5000 signal events at mu = 800 GeV in CompHep and slightly 
relaxing the jet transverse momentum cut, {Ptj > 50 GeV), 40000 background events in MadGraph were generated. 
The D quarks in signal events were made to decay in CompHep into SM particles. The final state particles for both 
signal and background events were fed into PYTHIA version 6.218 |21J for initial and final state radiation, as well 
as, hadronization using the CompHep-Pythia and MadGraph-Pythia interfaces provided by Athena (the ATLAS offline 
software framework) v9.0.3. To incorporate the detector effects, all event samples were processed through the ATLAS 
fast simulation tool, ATLFAST 1221 . and the final analysis has been done using physics objects that it produced. 

It must be noted that ATLFAST uses a parameterization for electrons, muons and hadrons without the detailed simulation 
of showers in the calorimeters. There is also a separate parameterization on the resolution for muons and electron tracks 
for the inner detector efficiencies. Minimum transverse energy of electromagnetic and hadronic clusters to be considered 
as electron or jet showers are Et > 5 GeV and Et > 10 GeV, respectively. Electromagnetic and hadronic cells in 
ATLFAST have the same granularity : A77 x A(p ~ 0.1 x 0.1. The electron isolation criteria requires a minimum distance 
AR = (Ayy)-^ x (Ac/))^ > 0.4 from other clusters and a maximum transverse energy deposition, Et < 10 GeV in 
outer cells accompanying the electron candidate. These outer cells are to be in a cone of radius AR = 0.2 along the 
direction of emission. The jets are reconstructed using the cone algorithm with the AR — 0.4 cone size. The smearing of 
particle clusters and jets is tuned to what is expected for the performance of the ATLAS detector from full simulation and 
reconstruction using the GEANT package 1231 . 

For this initial feasibility study, where the aim is to reconstruct the invariant mass of both D quarks, only e and /i decays 
of Z bosons are considered. Although the effective cross section becomes small compared to other decay channels, the 
benefits of this selection for a clean signal and for correct invariant mass reconstruction are indisputable. Other studies 
involving the invisible decays of Z boson and leptonic decays of W boson are in preparation. For the initial state particles, 
gg, uu and dd sub-channels were studied separately, and for final state particles only light quark jets were considered, all 
using CTEQ6L1 as PDF (TQ. For the case of AId — 800 GeV, the contributions from each sub-channel to the final cross 
section were about 50 %, 32% and 18% respectively. 

4.1 4/ 2jet channel, both Z fi^i case 

To reconstruct the two Z bosons, four isolated muons were required. The efficiency of finding 4 isolated muons is roughly 
50%. Since it is not known which two muons were originating from the same Z boson, all three possible combinations 
(assuming the lepton charge is not measured) were considered. The reconstructed invariant mass of the muon pairs was 
required to be close to the Z boson mass: within a window of 20 GeV around the central value of 90 GeV. If all pair 
combinations happened to be within the mass window, the one with the invariant mass closest to the measured value of Z 
boson mass was taken. The efficiency of reconstructing two Z bosons was about 90%. High transverse momentum cuts 
as in equation set below were applied to all leptons and jets to distinguish the signal events from background. One should 
note that, all the muons in such an event would satisfy the currently set trigger conditions for leptons Pr.iepton > 20 GeV 
II19I . The momentum distributions for all muons in an event and the imposed cuts for both signal and background can be 
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Figure 4: The transverse momentum cuts for muons (upper set) and jets (lower set). The plots for signal at M o = 800 GeV 
and the SM backgrounds are shown with arrows pointing at the cut values. The highest jet Pt for signal events peaks 
around 300 GeV, whereas for background no such peak is observed. 




Figure 5: The invariant mass plots for signal at Mb — 800 GeV and background events in the Z 
from two reconstructed Z bosons and two highest Pt jets. 
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Figure 6: The D quark invariant mass reconstruction in Z 
background. The signal was generated for AId = 800 GeV. 
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seen in Fig.0]upper set where the vertical arrow points at the imposed cut value. The percentage of muons surviving the 
transverse momentum cut in equation (|5} is about 60%. The events with at least two jets with transverse momentum greater 
than 100 GeV were kept in order to fully reconstruct the invariant mass of D quarks. In the lower set of Fig.|3 the solid 
line shows the momenta of the two most energetic jets in an event, whereas the dashed curve is for all other jets in the same 
event, again the vertical arrows pointing at the cut values. The two most energetic jets and the two previously reconstructed 
Z bosons were combined to reconstruct the D quark pair which was the goal of this study. However, the correct association 
of the two most energetic jets to the two reconstructed Z bosons is not known and involves combinatorics. Since the mass 
of the D quark is not known a priori, the "mass window around the central value" method that was used for Z boson 
reconstruction cannot be applied. The effect of wrong jet-Z association has the impact of enlarging the tails for the signal 
invariant mass distribution. This problem was partially solved by selecting the combination with the smallest absolute 
value of the difference between the two reconstructed D quark masses. Reconstructed invariant mass histograms are given 
separately in Fig.|5]for both signal and background cases. Therefore, the list of all the analysis level cuts becomes: 



Pt,^ > 40 GeV, (5) 
Mz = 90 ± 20 GeV , 

Njet > 2, 

Prjet > 100 GeV . 

The selection cut efficiencies for both signal and background events are given in Table|2]where each efficiency is calculated 
relative to the previous one, and the last column contains the combined efficiency value. 



Table 2: The individual selection cut efficiencies in percent for the both Z ^ case, for signal and background 



channel 




Mz 






Pr^^et 


^co7nbi7ied 


Signal 


48 


91 


59 


100 


95 


25 


Background 


34 


96 


16 


96 


12 


0.6 



4.2 4/ 2jet channel, both Z ^ ee case 

The cuts for this channel are same as for the channel above (Eq. (|5}), except the cut on the transverse momentum of 
electron, PT.e>15 GeV has been used. The existence of at least one electron originating from the Z boson decay which 
satisfies the isolated lepton trigger condition has also been checked for all 4 electron events. The method for reconstructing 
the Z boson and eventually the D quark remains the same as in the previous case. The selection cut efficiencies for this 
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Figure 7: The D quark invariant mass reconstruction in ZZ ee^ii case for all signal sub-channels as compared to SM 
background. The signal was generated for Mjj = 800 GeV. 




case are shown in Table |3l for both signal and background. The electron cuts leave slightly higher surviving events for 
both signal and background as compared to the muons. This difference and the different cut values for lepton transverse 
momenta can be attributed to the over optimistic electron reconstruction in the ATLAS software. The impact of this issue, 
which was under investigation at the time of this writing, was estimated by comparing the fits to the reconstructed Z 
bosons from both electrons and muons and was found to be about 5%. The invariant mass spectra is given in Fig.|6lfor a 
bin width of 40 GeV, showing the expected number of signal events being larger than the background ones. 



Table 3: The individual selection cut efficiencies in percent for the both Z ^ ee case for both signal and background. 



channel 




Mz 


Pt.c 




Prjet 


^combined 


Signal 


40 


99 


87 


100 


95 


33 


Background 


37 


98 


83 


94 


7 


2.0 



4.3 4/ 2jet channel, one Z ^ ee and one Z ji^i case 

This case is based on two isolated electrons, two isolated muons and two jets. There is no ambiguity in the lepton selection 
for Z invariant mass reconstruction thus a simpler reconstruction algorithm suffices. Both muons and at least one electron 
satisfy the trigger condition. The event selection cuts are summarized in the equation set|6l 



iV^ - 2, iVe==2, (6) 
Pt,m > 40 GeV, 
PT,e > 15 GeV , 

Njet > 2 

PT,jet > 100 GeV, 
Mz = 90 ± 20 GeV , 

The selection cut efficiencies are given in TableH Since the branching ratio is higher by a factor of two, compared to the 
first and second cases, this case yields more signal events and dominates the results. The reconstructed invariant mass for 
the signal and the SM background is given in Fig.0showing that the expected number of signal events is higher than for 
the background, in the region of the peak. 



5 Results 



All three above mentioned leptonic reconstruction cases were also considered for systematic studies at other D quark 
mass values;A/D = 600, 1000, 1200 GeV (The details of this analysis can be found in L20J ). For values of D quark mass 
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Table 4: The individual selection cut efficiencies for one Z ee and one Z ^ sub-case. The subscript I represents 
both electron and muon cases. 



channel 


Ni 


Mz 


Pt,i 




PT,jet 


^combined 


Signal 


44 


94 


71 


100 


93 


28 


Background 


35 


97 


34 


95 


10 


1.1 



larger than 900 GeV, in order to increase the expected statistical significance of signal identification, the cuts on the basic 
kinematic variables were modified as: 

Pt,m > 50 GeV, (7) 
Pr.e > 20 GeV , 
PTjet > 120 GeV . 

Using the convention of defining a running accelerator year as 1 x lO^seconds, one LHC year at the full design luminosity 
corresponds to 100 fh^^. For one such year worth of data, all the signal events are summed and compared to all SM 
background events as shown in Fig.|8l It is evident that for the lowest of the considered masses, the studied channel gives 
an easy detection possibility, whereas for the highest mass case (M£i=1200 GeV) the signal to background ratio is of the 
order of unity. If Nature has assigned a high mass to the D quark, other possible detection channels would be either to tag 
one Z via its leptonic decays and consider the neutrino decays of the other, or would involve hadronic decays of at least 
one Z and methods to disentangle the jet association. The detailed study of these modes as well as the charged current 
decay modes is deferred to future work. 

For each D quark mass value that was considered, a Gaussian is fitted to the invariant mass distribution around the D 
signal peak and a polynomial to the background invariant mass distribution. The number of accepted signal {S) and 
background {B) events are integrated using the fitted functions in a mass window whose width is equal to 2cr around the 
central value of the fitted Gaussian. The significance is then calculated at each mass value as S/^/B, using the number 
of integrated events in the respective mass windows. The expected signal significance for three years of nominal LHC 
luminosity running is shown in Fig. |9]left hand side. The shaded band in the same plot represents the statistical errors 
originating from the fact that for each signal mass value, a finite number of Monte Carlo events was generated at the 
start of the analysis and the surviving events were selected from this event pool. The statistical errors were calculated 
using binomial distribution and propagated to all significance and luminosity calculations. Given the small value of event 
yield, the Poisson probability distribution is more appropriate to set the observation and discovery limits 1241 . Following 
the PDG, the Gaussian probabilities for observing 2, 3 and 5 standard deviations when no signal is expected are defined 
as 4.55%, 0.27% and 5.7 x 1Q^^% respectively. Since the number of observed events in Poisson distribution (and also in 
experiment) is an integer, the values approaching their Gaussian counterparts from the lower side are taken. Consequently, 
this procedure always yields more pessimistic luminosity values. For example at M^i = 800 GeV, the probability utilized 
for 3(7 limit is 1.34 x 10~'^, about half of its Gaussian counterpart. Therefore, the reach of ATLAS to either exclude 
the existence of or discover the D quark using the studied dilepton channel is given in Table |5l for different D quark 
mass values. We observe that for A/d = 600 GeV, ATLAS could observe the D quark with a significance more than 3 
sigma before the end of the first year low luminosity running (10 /&^^/year) whereas to claim discovery with 5 sigma 
significance, it would need about 2 years of running time at low luminosity. For Mu — 1000 GeV, about two years of 
high luminosity running is necessary for a 3 sigma signal observation claim. The graphical representation is also given in 
Fig|9lright hand side. 

For these results, a possible source of systematic error is the selection of the QCD scale: although in a complete com- 
putation using all possible diagrams including loops this selection becomes irrelevant, this work relies on tree level cal- 
culations. An increase in QCD scale from the presently used value of Z boson mass up to D quark mass would mean 
a decrease in the total cross sections of about 65%. Such a change would also affect the results presented in this work, 
taking the 3 sigma signal observation limit mass from 1000 GeV down to 840 GeV for two years of high luminosity 
running. It should also be noted that the detector related backgrounds and the pile-up effects which could be important at 
high luminosity LHC, are not considered in this work. 



6 Conclusions 

This study shows that, for a range of D quark mass from 600 to 1000 GeV, ATLAS has a strong potential to find new 
physics related to GUT. After three years of design luminosity running, the 5 sigma discovery reach for D quark mass 
is about 920 GeV . The impact of the assumptions about D quark mass hierarchy and inter-family mixing is minimal: If 
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Figure 8: Combined results for possible signal observation at Mu ~ 600, 800, 1000, 1200 GeV . The reconstructed D 
quark mass and the relevant SM background are plotted for a luminosity of 100 fb^^ which corresponds to one year of 
nominal LHC operation. The dark Une shows the signal and background added, the dashed Une is for signal only and the 
light line shows the SM background. 
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Table 5: The required integrated luminosity in fb ^ to discover D quark as a function of its mass is shown. Expected 
signal and background event number are also given for one year of high luminosity running. 



D quark mass (GeV) 


600 


800 


1000 


1200 


signal & background events / year 


38.2 3.0 


9.6 1.3 


2.0 0.4 


0.6 0.2 


Luminosity for 2 a signal (fb^^) 


2.4 


18.4 


86 


373 


Luminosity for 3 a signal (fb^^) 


7.3 


36.7 


215 


870 


Luminosity for 5 a signal (fb^^) 


19.4 


91.7 


559 


2480 



the S quark is the lightest or the D quark mixes mostly to the second family, all the results and conclusions stay as they 
are. If the third family is involved, either B quark being the lightest or through large mixings to the third quark family, 
the results would also remain valid, provided no distinction of third family quarks is imposed. However, if the b quarks 
are to be identified to disentangle the D quark mixing, the h — jet tagging efficiencies should be convoluted with the 
presented results which would reduce the number of expected signal and background events by at least 50%. To enlarge 
the experimental reach window for the higher D quark mass values, the inclusion of other channels (remaining lines of 
the Table [0 is also envisaged. Furthermore, the inclusion of additional gauge bosons predicted by the Eq group could 
enhance the signal in the s channel if they have suitable masses. These studies are in preparation. 
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Figure 9: On the left: the expected statistical significance after 3 years of running at nominal LHC luminosity assuming 
Gaussian statistics. The vertical line shows the limit at which the event yield drops below 10 events. On the right; the 
integrated luminosities for 3 sigma observation and 5 sigma discovery cases as a function of D quark mass. The bands 
represent statistical uncertainties originating from finite MC sample size. 
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